function my_Laplacian()

image = imread("building.tif");
image_double = im2double(image);

%测试滤波
%Laplacian_mask = [1/16,1/8,1/16;1/8,1/4,1/8;1/16,1/8,1/16];
%测试锐化
%Laplacian_mask = [0,-1,0;-1,4,-1;0,-1,0];
Laplacian_mask = [-1,-1,-1;-1,8,-1;-1,-1,-1];
conv_result = my_conv(image,Laplacian_mask);
[row,col] = size(Laplacian_mask);
[m,n] = size(image);

for i=((row-1)/2+1):(m-(row-1)/2)
    for j=((col-1)/2+1):(n-(col-1)/2)
        image_double(i,j) = image_double(i,j) + conv_result(i-1,j-1);
    end
end

subplot(1,3,1);imshow(image),title("原始图像");
subplot(1,3,2);imshow(conv_result),title("轮廓");
subplot(1,3,3);imshow(image_double),title("变换后图像");

end